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Motivated by recent experiments of successfully carving out stable carbon atomic chains from graphene 
[Meyer et al Nature (London) 454, 319 (2008); Jin et al Phys. Rev. Lett. 102, 205501 (2009)] we inves- 
tigate a device structure of a carbon chain connecting two zigzag graphene nanoribbons with highly tunable 
spin-dependent transport properties. Our calculation based on the non-equilibrium Green's function approach 
combined with the density functional theory shows that the transport behavior is sensitive to the spin configura- 
tion of the leads and the contact position in the gap. A connection in the middle gives an overall good coupling 
except for around the Fermi energy where the leads with anti-parallel spins create a small transport gap while 
the leads with parallel spins give a finite density of states and induce an even-odd oscillation in conductance 
in terms of the number of atoms in the carbon chain. On the other hand, a connection at the edge shows a 
transport behavior associated with the spin-polarized edge states, presenting sharp pure a-spin and /3-spin peaks 
beside the Fermi energy in the transmission function. This makes it possible to realize on-chip interconnects or 
spintronic devices by simply changing the spin configuration of the leads and the position of the connection. 



I. INTRODUCTION 

Carbon based nanostructures, especially, quasi- ID struc- 
tures like carbon nanotubes (CNTs) and, recently, graphene 
nanoribbons (GNRs), are playing more and more important 
role in the development of nanoelectronics,^ possibly lead- 
ing to an era of carbon-based electronics. The practical ap- 
plication of CNTs is, however, limited by the challenge in 
controlling experimentally their diameter and chirality which 
control whether they are metallic or semiconducting. Differ- 
ent from CNTs, the electronic properties of GNRs are de- 
termined by their edge geometry and width and have shown 
promise for future generation of transistor.^"^^ Zigzag-edged 
GNRs (ZGNRs) are particularly intriguing because of their 
spin-polarized edge states which are localized around the two 
edges. The resulting spin-dependent transport may make them 
promising candidates for applications in spintronics. In nar- 
row ZGNRs with n zigzag carbon chains (denoted by n- 
ZGNRs, n ^ 32) the anti-parallel spin configuration is the 
ground state, which is slightly more favorable in energy than 
the parallel-spin one, while in wider n-ZGNRs (n ^ 32) the 
two spin configurations are both possible to exist due to the 
negligible interaction between the two edge states. Even in 
narrow ZGNRs, the small energy difference makes it possible 
to change the ground state from anti-parallel spins to parallel 
spins by appling a magnetic or electric field. 

Recently, an interesting progress related to graphene is the 
successful fabrication of free-standing linear carbon atomic 
chains carved out from a graphene sheet by high-energy elec- 
tron beam,^^"^"^ which is found stable and is connected by sp^ 
bonding. Unlike CNTs and GNRs, a carbon atomic chain has 
no chirality and width. Therefore, it provides an ideal trans- 
port channel for molecular devices. 

Experimentally, a carbon chain made in this way can be 
used as an on-chip device with the advantage of the perfect 
sp^ connection to the leads already set. This is in striking 
contrast to the situation in conventional molecular electron- 



ics using metal electrodes where a well-defined molecule-lead 
contact with a good reproducebility is a big challenge. For 
the miniatrization of the whole device, the two graphene leads 
can be cut to form GNRs. 

Theoretically, the device structure of a carbon chain con- 
necting two ZGNR leads is particularly interesting since it 
combines the simple transport channel with the rich electronic 
properties of ZGNRs via the perfect sp^ contact. Its transport 
properties can be artificially tailored to realize different func- 
tionalities, as shown in this work. First, the tunable spin state 
of the leads can be used to control the density of states at the 
Fermi energy, giving either a large or zero equilibrium con- 
ductance (switch). Second, the atomic thin bridge can be used 
to explore the locally spin-polarized edge states, implying that 
the transport behavior will be sensitive to the position of the 
bridge. This may be used to realize selectively the functional- 
ity of either on-chip interconnects or spintronics. 

So far, the spin-dependent transport phenomena in this kind 
of junctions has not been explored to the best of our knowl- 
edge, though the conductance of carbon chains connected to 
different types of electrodes have recently been studied by 
several research groups but how it behaves is still an 
open problem depending on the nature of the electrode used. 
Among these work, a calculation^^ of spin-denpendent trans- 
port was reported for a junction with very narrow ZGNR leads 
connected via a spin-polarized 5 -member ring which is, how- 
ever, unlikely to be carved out directly from a graphene sheet 
and is probably unstable due to the unpaired Pz electron with- 
out doping. Very recently Shen et al reported a theoretical 
calculation for such junctions with very narrow ZGNR leads 
but ignoring totally the spin freedom*2S An interesting find- 
ing from their calculation is that the equilibrium conductance 
shows an even-odd oscillation with regard to the number of 
atoms in the carbon chain^^: A chain with an odd number of 
atoms will have a larger conductance than one with an even 
number of atoms. However, if one takes the spin freedom into 
account, this even-odd behavior will totally disappear for the 
ground state since the anti-parallel spins will create a band gap 
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in the leads. 

In this work, we investigate the spin-dependent elec- 
tron transport in ZGNR-carbon chain junctions by perform- 
ing first-principles calculations based on the non-equilibrium 
Green's function (NEGF) approach combined with the den- 
sity functional theory (DFT)j22i22 We study systematically the 
effects of different spin configurations of the leads, the posi- 
tion of the carbon chain connection in the gap, as well as the 
width of the leads. Our calculation shows that the position of 
the bridge connection combined with the spin freedom of the 
leads dominates the transport behavior. 

The present work indicates that different functionalities can 
be achieved by changing the contact position and the spin 
state of the leads: on-chip interconnects can be realized with a 
bridge connected in the middle and a parallel spin configura- 
tion of the leads, which give a large equilibrium conductance, 
while a device for spintronics may be realized with a bridge 
connected at the edge, which gives a high spin-polarization 
ratio. 
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FIG. 1: Optimized atomic structures of the (4-ZGNR)-C7-(4-ZGNR) 
and (6-ZGNR)-C7-(6-ZGNR) junctions with different positions of 
connection: (a), (d) in the near middle (Ml -bridge), (b), (e) at the 
edge (E-bridge), and (c) in the right middle (M2-bridge). 



II. COMPUTATIONAL DETAILS 

We consider two widths for the ZGNR leads, consisting of 
4 and 6 zigzag carbon chains, respectively (see Fig. [T]). A 
carbon chain with 7 atoms is studied in details and a chain 
with 6 atoms is also calculated for showing the even-odd be- 
havior in the equilibrium conductance. For the junction with 
the 4-ZGNR leads, the carbon chain is connected at two po- 
sitions, one in the middle (Ml -bridge) and the other at the 
edge (E-bridge). For the junction with the 6-ZGNR leads, 
three connection positions are considered: in the right mid- 
dle (M2-bridge), in the near middle (Ml -bridge), and at the 
edge (E-bridge), respectively, (see Fig.[T]). 

Each edge of the two ZGNR leads has two possible spin 
states, a- or /3-spin. In this work, we consider two spin con- 
figurations for the leads. In one configuration, both the left 
and right leads are in anti-parallel spin state with the top edge 
being a- spin and the bottom edge being /3-spin (labeled by 
(a/3, a/3)), which is the ground state of the narrow ZGNR 
lead. In the other, the two leads are in parallel-spin state with 
all the four edges being a-spin (labeled by (aa, aa)), which is 
possible under a magnetic field and/or for wider ZGNR leads. 
Our calculation will show that the effect from the width of 
the lead is quite small and therefore the narrow ZGNR leads 
considered will also reflect the major behavior of wider ones. 

To investigate the electron transport through the 
ZGNR-carbon chain junctions we adopt the NEGF-DFT 
approach^^'^^ which combines the NEGF formula for trans- 
port with ab initio DFT calculation for electronic structure. In 
practice, the infinitely long 1-D system is divided into three 
parts: left lead, right lead, and device region containing the 
carbon chain plus enough ZGNR layers to accommodate the 
carbon chain-ZGNR interaction. The self-consistent Kohn- 
Sham Hamiltonian of the device region and the self-energies 
of two semi-infinite ZGNR leads are used to construct a 
single-particle Green's function from which the transmission 
coefficient at any energy is calculated. The conductance G 
then follows from a Landauer-type relation. The computa- 
tional techniques have been described in details previously. 
Briefly, for the DFT electronic structure calculation, we use a 
numerical basis set to expand the wave function^ ^: A double 
zeta plus polarization basis set (DZP) is adopted for all atomic 
species. The local density approximation (LDA)^^ is used 
for the electron exchange and correlation and the optimized 
TrouUier-Martins pseudopotentials^^ are used for the atomic 
cores. The atomic structure of the junctions including 
the carbon chain-ZGNR separation are fully optimized by 
minimizing the atomic forces on the atoms to be smaller than 
0.02 eV/A. 



III. RESULT AND DISCUSSION 

A. Effects of spin and position of bridge connection 

The calculated transmission functions for the (4-ZGNR)- 
C7-(4-ZGNR) and (6-ZGNR)-C7-(6-ZGNR) junctions are 
plotted in Figs. |2]and[3l respectively. The first thing to note is 
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FIG. 2: Transmission functions of the (4-ZGNR)-C7-(4-ZGNR) 
junction. The first panel: Ml-bridge connection in (a) (a[3, a[3) 
and (b) (aa, aa) spin configuration. The second panel: E-bridge 
connection in (c) (a/3, a/3) and (d) (aa, aa) spin configuration. 
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FIG. 3: Transmission functions of the (6-ZGNR)-C7-(6-ZGNR) 
junction. The first panel: M2-bridge connection in (a) (a/3, a/3) 
and (b) (aa, aa) spin configuration. The second panel: Ml-bridge 
connection in (c) (a/3, a/3) and (d) (aa, aa) spin configuration. The 
third panel: E-bridge connection in (e) (a/3, a/3) and (f) (aa, aa) 
spin configuration. 



that the transmission function depends very significantly on 
the spin configuration of the leads and the position of the 
carbon chain connection in the gap. When the connection 
is positioned around the middle range of the gap (M2- and 
Ml-bridge) the transmission function is overall broad in the 



energy window [-1, 1] eV, indicating a strong coupling be- 
tween the carbon chain and the ZGNR leads. For the anti- 
parallel spin configuration (a/3, a/3) there is a small transport 
gap around the Fermi energy (Figs. |2](a) and [3] (a), (c)) while 
for the parallel spin confugration (aa, aa) there are sharp res- 
onance peaks there (Figs. |2](b) and[3](b), (d)). The transport 
gap in the case of (a/3, a/3) is due to the band gap created by 
the anti-parallel spins in the leads, as shown in Figs. H] (a) and 
(c). Note that this band gap decreases slowly with the increas- 
ing width of the ribbon and so is the transport gap: 0.37 and 
0.35 eV for the 4-ZGNR and 6-ZGNR junctions, respectively. 
In the case of (aa, aa), the parallel spin state causes a band 
crossing at the Fermi energy as shown in Figs.|4](b) and (d), 
providing an finite DOS around the Fermi energy in the leads. 
The coupling of this finite DOS to the electronic states in the 
carbon chain gives rise to the sharp resonance peaks around 
the Fermi energy (see futher discussion later). 

On the other hand, when the carbon chain connection is 
positioned at the edge of the gap (E-bridge) the transmission 
function becomes very sharp peaks in the energy window [-1, 
1] eV (Figs. [2](c), (d), and[3](e), (f) ), indicating an over- 
all weak chain-lead coupling. Furthermore, around the Fermi 
energy the a- and /3-spin peaks are now separated from each 
other, making the transport at the Fermi energy is also blocked 
for the (aa, aa) spin configuration (Figs. |2](d) and[3](f) ). 

In the case of (aa, aa) the splitting between the a- and 
/3-spin peaks around the Fermi energy varies with the position 
of the bridge connection. For the wider 6-ZGNR junction it 
becomes more and more larger as the bridge is closer to the 
edge, from M2- to Ml- and then to E-bridge (see Figs. [3] (b) 
vs (d) vs (f)), because of the increasing effect from the spin- 
polarized edge states. 

As for the effect from the width of the lead, one can see that 
for the same Ml-bridge connection the overall effect in the 
transmission function is not significant and does not change 
the qualitative result when the width of the leads is increased 
from 4 to 6 zigzag chains: Fig. |2] (a) vs Fig. [3] (c) for (a/3, 
a/3), and Fig.[2](b) vs Fig.[3](d) for (aa, aa). Therefore, the 
results from the present calculation can be expected to be still 
qualitatively valid for junctions with wider ZGNR leads. 



B. Analyses 

Physically, the complicated spin-dependent transport be- 
havior and the very significant effect from the contact position 
is related to the spin-polarized edge states and their spatial 
distribution in the leads. To understand the results and pro- 
vide an insight into the physics underlying, we study the spin- 
polarized projected density of states (PDOS) for the leads and 
the local density of states (LDOS) for the characteristic peaks 
in the transmission function of the junctions. 

In Fig. [5] we show the spin-polarized PDOS for the edge and 
bulk regions of an infinite long 6-ZGNR. In the bulk region of 
the ribbon (Figs. [5] (a) and (b) ) the bulk states provide a finite 
DOS in the whole energy window except for that around the 
Fermi energy where a gap appears for (a, /3) because of the 
band gap shown in Fig. H] (c), while a finite DOS still exists 
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FIG. 4: The first panel: Band structure of a 4-ZGNR in (a) (a, /3) 
and (b) (a, a) spin configuration. The second panel: Band structure 
of a 6-ZGNR in (c) (a, /3) and (d) (a, a) spin configuration. 
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FIG. 5: The first panel: PDOS projected on the bulk region of a 6- 
ZGNR for (a) (a, 13) and (b) (a, a) spin configuration. The second 
panel: PDOS projected on the edge region of a 6-ZGNR for (c) (a, 
13) and (d) (a, a) spin configuration. The third panel: PDOS pro- 
jected on the C7 chain in the (6-ZGNR)-C7-(6-ZGNR) junction with 
(e) M2-bridge and (f) Ml -bridge connection in the {aa, aa) spin 
configuration. 

for (a, a) due to the band crossing shown in Fig.|4](d). When 
the bridge is positioned in the middle region of the gap the fi- 
nite DOS apart from the Fermi energy couples with the states 
in the bridge, leading to the very broad transmission function 
shown in Figs. [3] (a) - (d). To show this more clearly, we plot 



(a) 




FIG. 6: (a) and (b): LDOS of the a- spin component around -0.6 eV 
in Fig|3](a) and (c), respectively. 




FIG. 7: (a), (b), (c): LDOS of the /3-spin peaks around the Fermi 
energy (-0.04, 0.06, and 0.17 eV, respectively) in Fig|3](b), (d), and 
(f), respectively. 

in Fig. [6] the LDOS for the energy around -0.6 eV in Figs. [3] 
(a) and (c), respectively. It can be seen that in both cases the 
LDOS spreads out through the whole junction and distributes 
quite evenly in the leads, indicating that it is associated with 
the bulk states. In the case of (a, a) the finite DOS around 
the Fermi energy couples with the states in the bridge (their 
PDOS is plotted in Figs. [5] (e) and if)) and gives rise to the 
resonance peaks around the Fermi energy in the transmission 
function (Figs. [3] (b) and (d)). To show the nature of these res- 
onance peaks we plot their LDOS for the /3-spin component in 
Figs. [7] (a) and (b) for the M2- and Ml -bridge connections, re- 
spectively. One can see that when the connection is in the right 
middle (M2-bridge) the LDOS has large contribution from the 
bridge region and is distributed evenly throughout the leads, 
indicating that it is a result of the coupling between the states 
in the bridge and the bulk states in the leads. When the bridge 
is closer to the edge (Ml -bridge) now the LDOS has also large 
contribution from the edge states within the scattering region 
and causes the splitting between the a- and /3-spin compo- 
nents. 

On the other hand, for the zigzag chain of the edge (Figs. [5] 
(c) and (d) ) the PDOS reflects essentially the edge state which 
has two sharp peaks beside the Fermi energy, i.e., the occupied 
a-spin and unoccupied /3-spin states, respectively, placing a 
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FIG. 8: Transmission function of the spin-unpolarized (4-ZGNR)- 
C7-(4-ZGNR) junction. 

near-zero DOS elsewhere. Consequently, when the bridge 
is positioned at the edge it mainly couples with these edge 
states, leading to the transmission function with only very 
sharp peaks corresponding to the edge states, as shown in 
Figs. (e) and (f). The LDOS of the /3-spin peak in Fig. [3] 
(f) plotted in Figs. [7] (c) shows that the transport is indeed 
along the two edges. In this case, even for the (aa, aa) spin 
configuration the transport around the Fermi energy is also 
nearly blocked. Note that unlike the transport gap created by 
the anti-parallel spins in (a/3, a/3) case, this transport gap be- 
tween the a- and /3-spin conponents actually increases from 
the 4-ZGNR to the 6-ZGNR junction (0.26 and 0.33eV, re- 
spectively). This is because the a — a spin interaction in the 
narrow 4-ZGNR lead induces extra dispersions of the edge 
states, as shown in Fig.|4](b) vs (d). 

The broad and very narrow peaks in the transmission func- 
tion reflect the overall coupling strength between the carbon 
chain and the ZGNR leads. It is interesting to note that the 
best coupling is given by neither the bridge at the edge (E- 
bridge) nor the bridge right in the middle (M2-bridge) but the 
one positioned in between (Ml -bridge). This is evident in the 
LDOS shown in Fig. [6] where the Ml -bridge gives a remark- 
ably larger LDOS around the bridge and contact region than 
the M2-bridge does. The result is a much larger transmission 
coefficient around -0.6 eV in Fig.[3](c) than that in Fig.[3](a). 



C. Even-odd behavior 

Finally, we would like to discuss the even-odd behavior 
in the equilibrium conductance found previously in a spin- 
unpolarized calculation.^^ This calculation shows that short 
carbon atomic chains with an odd number of atom will give 
a significantly larger equilibirum conductance than those with 
an even number of atom. According to the present calcula- 
tion, however, the spin freedom must be taken into account 
in reaching this conclusion. To have a direct comparison, 
we also perform spin-unpolarized calculation for the narrow 
4-ZGNR junction with the Ml -bridge and plot the result in 
Fig.[8l which is basically the same as obtained in Ref..^^ The 
very sharp peak at the Fermi energy is the origination of the 
even-odd oscillation in conductance, which was ascribed to 
the edge states. The present spin-polarized result in Fig. |2] 
shows that this peak around the Fermi energy only exists in the 
case of (aa, aa) with the Ml-bridge connection (see Fig.|2] 



FIG. 9: (a) Transmission function of the (6-ZGNR)-C6 -(6-ZGNR) 
junction in the (aa, aa) spin configuration and with the M2-bridge 
connection, (b) PDOS projected on the Ce chain of the junction. 

(b)) but now having a small splitting between the a- and /3- 
spin components. In the case of (a/3, a/3) the equilibrium 
conductance is zero simply because the band gap created by 
the anti-parallel spins in the leads (see Fig. [2] (a)). Addition- 
ally, for the E-bridge connection even in the case of (aa, aa) 
the equilibrium conductance is also nearly zero because of 
the large splitting between the a- and /3-spin components (see 
Fig.|2](d)). For the wider (6-ZGNR)-C7-(6-ZGNR) junction, 
the result is similar - the sharp peaks around the Fermi energy 
only appear in the case of (aa, aa) and only when the bridge 
is positioned around the middle region (see Fig. [5]). 

The present spin-polarized calcuation shows that, as can be 
seen in Figs. [7] (a) and (b), the large resonance peaks at the 
Fermi energy is not originated from the edge states but comes 
from the coupling between the finite DOS in the leads due to 
the band crossing and the states in the bridge whose PDOS is 
given in Figs. [5] (e) and (f) showing a large DOS around the 
Fermi energy for the C7 chain. Since in this case the coupling 
to the leads is determined by the states in the carbon chain, the 
resulting equilibrium conductance will be sensitive to its elec- 
tronic state which may be affected significantly by the number 
of atoms in the chain. As was found in Ref.,^^ an even or odd 
number of atoms in the carbon chain will result in different 
C-C bond-length distribution due to the Peils transition effect, 
which gives quite different electronic states. 

In order to make this issue more clear, we calculate the 
transimission function of a (6-ZGNR)-C6-(6-ZGNR) junction 
with the M2-bridge connection for the (aa, aa) spin config- 
uration. The result and the PDOS projected on the Ce chain 
are plotted in Figs.|9](a) and (b), respectively. The small DOS 
at Fermi energy for the Ce chain shows the the coupling be- 
tween the lead and the carbon chain is very weak, resulting a 
much smaller equilibrium conductance compared with the C7 
junction. The consequence is an even-odd oscillation in the 
equilibrium conductance in terms of the number of atoms in 
the carbon chain. 



IV. SUMMARY 

In summary, motivated by recent experiments of suc- 
cessfully carving out stable carbon atomic chains from 
graphene, we have studied the spin-dependent electron trans- 
port through carbon atomic chains connecting two zigzag 
graphene nanoribbons by using the non-equilirium Green's 
function approach combined with the density functional the- 
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ory calculation. The effects on the transport from different 
spin configurations of the leads, different positions of the 
bridge connection, and the width of the leads are investigated. 

It was found that the spin freedom combined with the po- 
sition of the bridge connection dominate the transport proper- 
ties. A bridge connection in the middle will give an overall 
good coupling and transparency because of the strong cou- 
pling between the bulk states in the leads and those in the 
atomic chain, except for the energies around the Fermi en- 
ergy where the lead with anti-parallel spins creates a transport 
gap while the lead with parallel spins give an finite density of 
states because of its band crossing. The coupling between this 
finite density of states and the states in the carbon chain gives 
rise to a large (for odd number of atoms) or a small (for an 
even number of atoms) equilibrium conductance for both the 
a- and /3-spin components, inducing an even-odd oscillation 
in the equilibrium conductance. On the other hand, a bridge at 
the edge leads to a transport behavior revealing the edge states, 
showing only sharp pure a- spin and /3-spin peaks beside the 



Fermi energy with a near zero equilibrium conductance. 

Our calculation reveals that a functional device for on-chip 
interconnects can be realized by a bridge connection in the 
middle, which gives a large equilibrium conductance (~ Go) 
when the leads are in parallel spin state. On the other hand, a 
functional spintronic device with a high spin polarization ratio 
may be realized by a connection at the edge with a small gate 
voltage shifting the Fermi energy to the energy of the pure a- 
or /3-spin peak. 
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